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Abstract. The conformal method for constructing initial data for Einstein's equations 
is presented in both the Hamiltonian and Lagrangian picture (extrinsic curvature decom- 
position and conformal thin sandwich formalism, respectively), and advantages duo to 
the recent introduction of a weight-function in the extrinsic curvature decomposition 
, are discussed. I then describe recent progress in numerical techniques to solve the re- 

suiting elliptic equations, and explore innovative approaches toward the construction of 
' astrophysically realistic initial data for binary black hole simulations. 
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1. Introduction 



■ Numerical methods play an important role for investigations into the properties 

of Einstein's equations. In particular, the late stages of inspiral and coalescence of 
binary compact objects like binary black holes are thought to be accessible only to 
numerical investigations. Knowledge of the full waveform of inspiraling binary black 
holes, including the highly nonlinear coalescence phase, will enhance sensitivity of 
gravitational wave detectors like LIGO or GEO600 through cross-correlation of 
the observed signal with the expected waveforms 20 . Comparison of the observed 
signals with the predictions of general relativity will test general relativity in the 
genuinely nonlinear regime. Besides the experimental urgency, the binary black hole 
problem is arguably the most fundamental dynamical problem in general relativity; 
however, it remains unsolved. 

Initial data forms the starting point for any evolution. For Einstein's equations, 
the most widely used method to construct initial data is the conformal method, 
pioneered by Lichnerowicz |30| and extended to a more general form by York and 
coworkers |49I34I52] . In two recent papers, York |^ and Pfeiffer & York com- 
pleted the conformal method: It is now available in a Lagrangian and in a Hamil- 
tonian picture (referred to as the conformal thin sandwich formalism and the ex- 
trinsic curvature decomposition, respectively), and both pictures completely agree. 
The transverse-tracefrcc part of the extrinsic curvature is now defined such that it 
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vanishes for any stationary spacetime. The method is now completely invariant to 
conformal transformations of the free data. 

The conformal method results in a set of coupled nonlinear three-dimensional el- 
liptic partial differential equations. Over the last few years, numerical techniques for 
solving these coupled elliptic equations were improved tremendously. Construction 
of binary black hole initial data is no longer limited by numerical capabilities, but 
by the incomplete understanding of the choice of free data and boundary conditions 
for the elliptic equations. 

Here, we present the Lagrangian and Hamiltonian pictures of the conformal 
method, including many details which may have been mentioned in passing in tech- 
nical papers, but were never presented in a coherent fashion. In the second part 
of this paper, we describe briefly numerical methods and then explore recent inno- 
vative approaches to the construction of astrophysically realistic binary black hole 
initial data. Throughout this paper, emphasis is placed on physical and numerical 
issues, rather than mathematical proofs. 

2. The initial value problem 

Using the standard 3+1 decomposition |ll52j of Einstein's equations, we foliate 
spacetime with spacelike t =const. hypersurfaces. Each such hypersurface has a 
future pointing unit-normal ti^, induced metric = '•''^g/ijy + nfj,n^ and extrinsic 
curvature if^^ = —■^Cngfj.i^- The spacetime metric can be written as 

ds^ = -N^dt + gij {dx' + P'dt) {dx^ + 13^ dt) , (2.1) 

where N and denote the lapse function and shift vector, respectively. N measures 
the proper separation between neighboring hypersurfaces along the surface normals 
and /3* determines how the coordinate labels move between hypersurfaces: Points 
along the integral curves of the time-vector = Nn^^ + f3^ (where f]^ — [0, 
have the same spatial coordinates x'. 

Einstein's equations decompose into evolution equations and constraint equa- 
tions for the quantities gij and Kij . The evolution equations determine how gij and 
Kij are related between neighboring hypersurfaces, 

^tg^J = -2NK,j + V,/3j + V^A (2.2) 
dtK,, = N {R,j - 2K,kK!y + KK,j - SnGS.j + inGg^.iS - p)) 

- V,V,iV + /3^-Vfci^„- + K,kV,l3^ + KkjV.p''. (2.3) 

Here, Vi and R are the covariant derivative and the scalar curvature (trace of the 
Ricci tensor) of gij, respectively, and K — Kijg^^ denotes the mean curvature. 
Furthermore, G stands for Newton's constant, p and Sij are matter density and 
stress tensor, respectively, and S = Sijg^^ denotes the trace of Sij. 

The constraint equations are conditions within each hypersurface alone, ensur- 
ing that the three-dimensional surface can be embedded into the four-dimensional 



February 7, 2008 3:11 WSPC/INSTRUCTION FILE Pfeiffer 



The initial value problem in numerical relativity 3 

spacetirae: 

R + K^ - K,jK''^ = IQiiGp, (2.4) 
V, {K^^ - g'^K) = SnGf, (2.5) 

with denoting the matter momentum density. Equation H2.4|l is called the Hamil- 
tonian constraint, and Eq. (|2.5(l is the momentum constraint. 

Cauchy initial data for Einstein's equations consists of {gij,K'^^) on one hyper- 
surface satisfying the constraint equations l|2.4|l and H2.5|l . After choosing lapse and 
shift (which are arbitrary and merely choose a specific coordinate system), Eqs. H2.2|l 
and H2.3|l determine {gij,K'^^) at later times. Analytically, the constraints equations 
are preserved under the evolution. In practice, however, during numerical evolution 
of Eqs. (|2.2|l and (|2.3|l or any other formulation of Einstein's equations, many prob- 
lems arise. 

The constraints 1)2. 4|l and (|2.5|l restrict four of the twelve degrees of freedom of 
{gij, K''^). As these equations are not of any standard mathematical form, it is not 
obvious which four degrees of freedom are restricted. Hence, finding any solutions 
is not trivial, and it is even harder to construct specific solutions that represent 
certain astrophysically relevant situations like a binary black hole. 

2.1. Preliminaries 

Both Hamiltonian and Lagrangian viewpoints use a conformal transformation on 
the spatial metric, 

9i3 = ^^~9ij (2-6) 

with strictly positive conformal factor ip. gij is referred to as the conformal metric. 
From 1)2. 6|l it follows that the Christoffcl symbols of the physical and conformal 
metrics are related by 

F}, - t), + 2^-1 (^ja^V + Sid.i; - ~g,kfdi^) , (2.7) 

which in turn implies that the scalar curvatures of gij and gij are related by 

R^'ip-'^R~8i>^^V^ip. (2.8) 

Equations were already known to Eisenhart ^19) . Furthermore, for any 

symmetric tracefree tensor S'*-', 

= i^'^^VjS'^ (2.9) 

where V is the covariant derivative of gij. Lichnerowicz [301 used Eqs. H2.6|l to H2.9|l 
to treat the initial value problem on maximal slices, K = 0. For non-maximal slices, 
the extrinsic curvature is split into trace and tracefree parts |34 |. 

K'^ ^ A'^ + -g'^K. (2.10) 
3 
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With (|2.8() and (|2.10|l . the Hamiltonian constraint 12. 4() becomes 

- -i^R - —ip^K^ + -iJ^A.jA'^ + 27rGV'V 0, (2.11) 
8 12 8 

a quasi-hnear Laplace equation for -0. Local uniqueness proofs of equations like 
H2.11|l often linearize around an (assumed) solution, and apply the maximum prin- 
ciple. However, the signs of the last two terms of (|2.11l) are such that the maximum 
principle cannot be applied and consequently, it is not immediately guaranteed that 
Eq. H2.11|l has locally unique solutions. The term proportional to AijA^^ will be dealt 
with later; for the matter terms we follow York and introduce conformally scaled 
source terms: 

/ = (2.12) 
p = ^-^p. (2.13) 

The scaling for f makes the momentum constraint below somewhat nicer; the 
scalings of p and are tied together such that the dominant energy condition 
preserves sign: 

- 9^jfj' = V--^' (p' - ~9^^fr) > 0. (2.14) 

With Eq. (|2.13|) . the matter term in l|2.11|) becomes 2ttGiP~^p with negative semi- 
definite linearization for p > 0. 

The decomposition of Kij into trace and tracefree part, Eq. H2.10|l . turns the 
momentum constraint H2.5() into 

VjA'^ -^V'K ^SirGf. (2.15) 

o 

The conformal transformation H2.6|l implies one additional conformal scaling 
relation. The longitudinal operator |16l5UI5l] 

{hVy^ = V'V^ + \/^V' - \g''VkV^, (2.16) 

satisfies [50) 

(i^vy^ =i>-^{lvy'. (2.17) 

Here {hVy^ is given by the same formula (|2.1()|l but with quantities associated with 
the conformal metric gij. (In fluid dynamics (LT^)*-' is twice the shear of the velocity 
field V^). In d spatial dimensions, the factor 2/3 in Eq. (|2.16|) is replaced by 2/d\ 
Eq. (|?T7|l holds for all d. 

2.2. Lagrangian picture — Conformal thin sandwich formalism 

The conformal thin sandwich formalism|53j deals with the conformal metric and its 
time derivative] as illustrated in Figure ^ we deal with two hypersurfaces separated 
by an infinitesimal St (explaining the name "thin sandwich" ) , and connected by 
lapse N and shift The mean curvature of each hypersurface is given by K and 
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tdt^{Nh + (3) 6t 




Fig. 1. Setup for conformal thin sandwich formahsm. 



K + dtKSt, respectively, and the metric is decomposed into conformal factor and 
conformal metric. This decomposition is synchronized between the two hypersur- 
faces by the requirement that the conformal metrics on both hypersurfaces have the 
same determinant to first order in St. The variation of the determinant of gij is 

S~9^~g~9''5~g,3^9~9''u,jSt, (2.18) 

so that Uij = dtcjij must be traceless. 

Besides the relationships indicated in Figure ^ the conformal thin sandwich 
formahsm rests on the nontrivial scahng behavior of the lapse function: 

N^il:^N. (2.19) 

Indications suggesting this scahng appear in a bewildering variety of contexts (see 
discussion in |5:-{I4()| 1. This scaling is also crucial in the present context as well, cf. 
Eq. ifT^ below. 

Substitution of Eq. H2.6(l into the evolution equation for the metric, Eq. H2.2|l . 
and splitting into trace and trace-free parts with respect to the physical inverse 
metric g''^ = tp~'^g'^^ results in 

dtln^p = -Ink + lVk(3'', (2.20) 



6 6 



^^u,j = -2NA,, + (L/3),,. (2.21) 



Ui-j = dtgij ~ ^gijg'^'dtgki- (2.22) 



Equation H2.21|l is the tracefree piece of dtgij, thus for Uij = ■ip'^Uij, 

1 
3 

We solve Eq. f^TT^i for A'^ , 

A^^ ^ ^ ((L/3)'^- - u^^) , (2.23) 
and rewrite with conformal quantities [using H2.19|l . (|2.17(l and u^^ = ^jj^^it^^]: 

A'i = ^-1" (^(L/3)*^' - u'i^ EE (2.24) 
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which defines the conformal tracefree extrinsic curvature A'-' . Equation (|2.24() shows 
that the formula for A^^ is form invariant under conformal transformations; this 
hinges on the scaling of N in Eq. (|2.19() . Substitution of Eq. (|2.24l) into the momen- 
tum constraint H2.15() and application of Eq. (|2.9|1 yields 

V, - V, {^-^n^') - ^^'V^^ = 87rGJ\ (2.25) 

whereas Eq. (|2.24(l modifies the Hamiltonian constraint (|2.11(l to 

- ivi? - ^V-'if' + \'^r^~A^,A'' - -2^G^-3 (2.26) 

O Iz <3 

Equations (|2.25l) and (|2.26l) constitute elliptic equations for /3* and "0. We can 
therefore construct a valid initial data set as follows: Choose the free data 

(5y-, uy-, K, N) (2.27) 

(and matter terms if applicable), solve Eqs. 12.25|) and (|2.26|) for and -0, and 
finally, assemble = -^^j and i^*^' = + \a'^K- 

We now comment on several issues related to the conformal thin sandwich for- 
malism. 



2.2.1. Fixing N via dtK 

In the free data Eq. (|2.27|) . ijij and = dtgij constitute a "variable & velocity pair" 
{q,q) in the spirit of Lagrangian mechanics, but the remaining free data does not. 
To improve this situation, we note that the trace of the evolution equation ()2.3|l . 
results in (see, e.g., {45 j ) 

dtK - (3^dkK = N{R + K"^ + 47rG'(5 - 3p)) - V^TV. (2.28) 

Elimination of R with the Hamiltonian constraint H2.4|l and rewriting this equation 
with conformal quantities results in 



V^A^ + 14V' In V'V.iV + N 

-h42Vj In V'V' In - AnGil^^iS+'ip) 



For given dtK, this constitutes an elliptic equation for N. Therefore, if we take 
dtK as the "free" quantity instead of N, then the free data for the conformal thin 
sandwich formalism becomes 

(g,,, u,j, K, dtK) (2.30) 

plus matter terms if applicable. These free data consist completely of (q, q) pairs as 
appropriate for the Lagrangian viewpoint. These free data are also useful in practice 
for computations of quasi-equilibrium initial data, for which dtK = is a natural 
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and simple choice, whereas it is not obvious which conformal lapse N should be 
used. We note that Eq. H2.29|l can be rewritten as 



ii? + l^^K' + + 27rGV/(p + 2S) 



(2.31) 



{dtK~/3''dkK) . 



This equation is somewhat shorter and computationally somewhat more convenient. 
Using the free data (|2.30l) . five coupled elliptic equations have to be solved, rather 
than four, namely and either or l(?3T)l . 



2.2.2. Invariance to conformal transformations of the free data 

Given free data l|2.27|) and a solution (■)/', of the conformal thin sandwich equa- 
tions, choose a function > 0, and define conformally rescaled free data by 

g,;^. M'^^ ^^-^-y, x' = K, N' = -i'-^N, (2.32) 

plus the scalings p' = '^^p, j' ' = \I>i"j ' for matter terms if applicable. These free 
data, together with conformal factor = ^'■0 and the shift /?' * = lead to the 
same physical initial data {gij , if *-' ) : 

9i3 = V-"* .9y- = ^J^9i] = gij, (2.33) 
^'y =^'-io_l. (^(L'/3')'^ -u'*^) =A'^. (2.34) 

Here, we used Eqs. (|2.6|l and (|2.24ll . and L' denotes the longitudinal operator of 5^, 
which, by Eq. (|TT7|) . satisfies (L'/3)'^ = ^'''(L/?)'^ . Adding the trace of the extrinsic 
curvature to Eq. (|2.34(l is trivial. 

Therefore, only the conformal equivalence class of gij is relevant for the physical 
solution. This is a very desirable property; we introduced gij as a conformal metric, 
so its overall scaling should not matter. Specification of dfK instead of N as part of 
the free data preserves this invariance, as the StX-equation, (|2.29(l . is derived from 
physical quantities in the first place. 

The extrinsic curvature decomposition introduced in the next section is also 
invariant under conformal transformations of the free data. We note that invariance 
under conformal transformations of the free data is not trivial; earlier variants of 
the constraint decompositions did not possess it, giving rise to ambiguities in the 
free data which influenced numerical investigations )45i] . 



2.2.3. Gauge degrees of freedom 

The physical initial data {gij, K'^^) has twelve degrees of freedom (matter just adds 
four additional degrees of freedom in p and J* which determine the four physical 
matter variables p and j*). It is constrained by four constraint equations, so there 
should be eight degrees of freedom in the freely specifiable data. However, even 
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taking into account that only the conformal equivalence class of gij is relevant and 
that Uij is traceless, the free data H2.27II or (12.30(1 consists of twelve quantities, not 
eight. To clarify this issue, consider the substitutions 

. ' ' (2.35) 

^ f}' + W\ 

The vector disappears from Eqs. H2.23|l - H2.26|l . therefore the substitution ((2.35(1 
will not change the physical initial data set {gij,K''^); it merely tilts the time-axis 
and changes the coordinate labels on the second hypersurface. Thus, Uij contains 
three gauge degrees of freedom associated with the shift. {W^ enters into Eq. ((2.29(1 
as an advection term, though, because dtK is the derivative along the time- vector) . 

The fourth "missing" degree of freedom is hidden in the lapse function N: One 
can construct every possible initial data set {gij , K^^ ) with any (non-pathologic) 
choice of N. This can be seen by going backward from the physical initial data 
{gij,K'^^) (satisfying the constraints) to the free data. Given {gij,K^^) and any N 
and /3' , set the free data ((2.27(1 by 

g^J=g^J, u'^ ^{Lpy^ -2NA'\ K = g,, (2.36) 

as well as the given TV. With these free data, ip = 1 and the given /?* will recon- 
struct the physical spacetime {gij,K^^) as can be seen from Eqs. ((2.6(1 and 112.24(1 . 
Therefore, -0 = 1 and the given /3* will solve the conformal thin sandwich equations 

and xrmt . 

The fact that we were free to choose /3* reflects again the gauge-symmetry il- 
lustrated in Eq. ((2.35(1 . but in addition, we showed that the choice of N does not 
restrict the set of "reachable" initial data sets. 

Note that the physical initial data contain further gauge freedom: Covariance 
under spatial transformations implies that g^j (and gij) contain three gauge degrees 
of freedom associated with the choice of coordinates. Furthermore, K can be inter- 
preted as time 49 , fixing the temporal gauge. Thus, an initial data set has only 
four physical degrees of freedom — in perturbed flat space they are simply the two 
polarizations of gravitational waves. 



2.2.4. Implications for an evolution of the initial data 

During the solution of the conformal thin sandwich equations, one finds a shift 
and a lapse N. If this gauge is used in a subsequent evolution of the initial data 
{gij,K'^^) then Eq. 1(2.22(1 implies that, initially ^ 

dtgij - \^9ijg''^dtgki = ip^Uij- (2.37) 

The freely specifiable piece Uij thus directly controls the tracefree part of the time- 
derivative of the metric. If we specified dtK as part of the free data, then, of course, 
this will be the initial time-derivative of the mean curvature. Finally, from ((2.20(1 . 
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we find 

dt\ni;^^{Vkj3'' ^NK) 

= lVkP''+P''dkln^P-li;^NK. (2.38) 
6 6 

Since we have been very successful so far with specification of time-derivatives 

{dtgij and dfK), one might be tempted to turn H2.38|l around and use it as the 

definition of K in terms of dt In V'. This idea certainly comes to mind when looking 

for quasi-equilibrium solutions, for which time-derivatives of as many quantities as 

possible should vanish. Pursuing this idea, we find from Eq. I|2.38|l 

K=-^(Vkf3'' -6 {dt -fS^dk) In • (2.39) 

Substituting Eq. (|2.39() into the momentum constraint Eq. H2.25|l . however, makes 
the principal part of Eq. H2.25|l proportional to 

dj&P' -d'dkP^, (2.40) 

which is non-invertible. Therefore the attempt to fix K via H2.39|l will fail. Indeed, 
during the construction quasi-equilibrium initial data sets of spherically symmetric 
spacetimes |15) . it was found that the conformal thin sandwich formalism with 
Uij = and dtK = Q (and appropriate boundary conditions) is so successful in 
picking out the time-like Killing vector that several different choices of K lead to 
solutions satisfying dt In?/' = 0. 

We conclude that, in contrast to the trace- free part H2.37|l . one can not easily 
control dt In ?/> by choices of the free data. One can only evaluate (|2.38|l after solving 
the conformal thin sandwich equations. 

2.3. Hamiltonian picture — Extrinsic curvature decomposition 

The second method to construct solutions of the constraint equations if based on 
a decomposition of the extrinsic curvature. Early variants of this approach 34 52] 
have been widely used for almost thirty years, but the final version was developed 
only very recently We will make use of the equations and results from section 
12.11 in particular, we use a conformal metric, gij — ip'^gij, and split the extrinsic 
curvature into trace and trace-free parts, K^^ = A^^ + l/3g^^K, cf. Eqs. H2.6|l and 

We start with a weighted transverse traceless decomposition of A^^ , 

A'^ = At^r^ + - {hvy' . (2.41) 
a 

Here, Al^j, is transverse, VjAl^j, — 0, and traceless, gijAl^j, ~ 0, and is a strictly 
positive and bounded function. Appearance of the weight function cr is a key point 
in the extrinsic curvature formulation; its inclusion is the major difference of |40| 
over the older variants. 
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Given a symmetric tracefree tensor like A^^ , the decomposition (|2.41f) is obtained 
by taking the divergence of Eq. H2.41|l . 

VjA'^ = Vj [a-\LVy^ . (2.42) 

The right hand side, Vj [cr^^(L .)*■'] , is a well-behaved elliptic operator in diver- 
gence form, so no problem should arise when solving H2.42(l for V^. Substitution of 
the solution back into 12.41|l yields Ajlrp. In the presence of boundaries, Eq. H2.42(l 
requires boundary conditions which will influence the solution and the decompo- 
sition H2.41(l . For closed manifolds, existence and uniqueness of the decomposition 
for the case tJ = 1 was shown in j51| . 
We now conformally scale the quantities on the right hand side of H2.41|l with 
the goal of rewriting the momentum constraint in conformal space. First, we set 

A!}j. = i;-^°Ai^rj.. (2.43) 

Equation H2.9|l ensures that Al^j, is transverse with respect to gij if and only if 
Alfrp is transverse with respect to the physical metric gij. Because of Eq. H2.17|l . 
and because L is the conformal Killing operator, the vector is not rescaled. The 
conformal scaling of the weight function is given by 

a = if'd. (2.44) 

The most immediate reason for this scaling is to allow Eq. H2.45|l below; several 
more reasons will be mentioned in the sequel. 

Using the scaling relations l|2.17|l . H2.43|l and H2.44|l . we can rewrite Eq. (|2.41|) 

as 



A'^ = '^-^^ ( A%T. ) - i^'^^A^ (2.45) 



where 



i'^' =i^y + a-i(LF)*^' (2.46) 



is a weighted transverse traceless decomposition in the conformal space. The scaling 
Ai = ip~^'^A'^^ was not postulated (as it had to be in the old variants), but follows 
from the other scalings. By virtue of the scaling of the weight function a, Eq. H2.44|l . 
the weighted transverse traceless decomposition thus commutes with the conformal 
transformation. This commutation of conformal transformation and weighted trans- 
verse traceless decomposition is precisely the new feature of the weighted decom- 
position. Without the weight-function, conformal transformations and transverse- 
traceless decomposition do not commute, leading to the two inequivalent old vari- 
ants, depending on which operation is performed first. 

Equations (j2.9|l and H2.45|l allow us to rewrite the momentum constraint H2.15|l 

as 

Vj (^^{LVy^^ ~ '^i^'^V'K = 8ttGj\ (2.47) 
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an elliptic equation for V^. The Hamiltonian constraint Eq. (|2.11|l reads 

1^, 1 
12 



- o-^^ - — V'^^^ + -i^^'^A.jA'^ = -2TrGip-^p, (2.48) 



with A^^ given by Eq. 12.46|l . Equation (|2.48l) is identical to Eq. I|2.26ll since in both 
formulations A^^ = ip~^^A^^ , however, the definitions of A^^ differ. 

Starting from the physical initial data {gij , K^^ ) , we have now rewritten the con- 
straints (|2.4|l and (|2.5() as elliptic equations H2.47|l and p.48|l . In order to construct 
a valid initial data set {gij, K'^-'), one first chooses the free data 

K, i^'y, a) (2.49) 

and matter terms if applicable, then solves Eqs. H2.47() and (|2.48() for and '0, and 
finally assembles the physical solution by Eqs. (|2.6() . H2.10|l . and H2.45|l . 



2.3.1. Remarks on the extrinsic curvature decomposition 

Similar to section 12.2.21 one can show that the physical initial data (gij , K'^^ ) is 
invariant to a conformal transformation of the free data. For 5* > 0, the relevant 
transformations are [cf. Eq. 



Ai^^^^'^'Ai^^, K'^K, a' = ^-'d, (2.50) 

plus the scalings p' — ^'^p, J'* — ^f^^J* for matter terms if applicable. The cal- 
culation is straightforward, the key-point being that the scaling of the weight- 
function H2.44|l synchronizes the conformal scaling of the transverse-traceless and 
longitudinal parts of the weighted transverse traceless decomposition. 

Because of the invariance to conformal scalings of the free data, gij supplies only 
five degrees of freedom, so that the free data Eq. I|2.49|l contains nine degrees of 
freedom. The weight a (or a) merely parametrizes the transverse traceless decom- 
position H2.41|l . For any choice of a, the decomposition 12.41|l can be performed, 
therefore with any choice of a, all initial data sets can be generated for appropriate 
choices of the free data. 

To construct a transverse traceless tensor Al^j, compatible with the metric gij, 
one decomposes a general symmetric tracefree tensor M*-' . Write 

M^' = + -(Lwy\ (2.51) 
a 

The divergence of this equation. 



a~'{hWy^ , (2.52) 



represents an elliptic equation for W\ Solving this equation, and substituting 
back into (|2.51() yields 

Ai^j, = - \{LWY^. (2.53) 
a 
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The formula for A^^ , Eq. (|2.46l) , now reads 

A'^ = M'^ + i [1{V - W)] , (2.54) 

which depends only on the difference — W\ On the other hand, subtraction of 
H2.52|l from the momentum constraint (|2.47l) yields 

Vj (^^[HV - W)]'^^ + Vj£'P^ - ^t/'^V'X = 8ttGj\ (2.55) 

which is an equation for the difference — W^. Thus one can combine the construc- 
tion of Al^j, from M^^ with the solution of the momentum constraint, as observed 
by Cantor HHj. Instead of solving (|2.52(l for and then 12.47|l for V\ one can 
directly solve l|?35l) for V - W\ 

In the presence of boundaries, solutions of elliptic equations like H2.55|l or H2.52(l 
will depend on boundary conditions. When constructing black hole initial data, in- 
ner boundaries are often present, and it is far from clear what boundary conditions 
to apply there. In Ref. 36 , for example, the situation is encountered that bound- 
ary conditions for the combined solution for V'^ — are known, but not for the 
individual solutions for and . 



2.3.2. Identification of a with the lapse N 



The extrinsic curvature formulation of the initial value problem as presented so far 
is perfectly adequate for the mathematical task of rewriting the constraints as well- 
defined equations. However, it is very natural to further identify the weight-function 
cr with the lapse function N , 



27V, 



2N. 



(2.56) 



One reason for this identification is that a and N have the same conformal scaling 
behavior, cf. Eqs. (|2.19() and H2.44|l . A second reason is that with this identification, 
the conformal thin sandwich equations become equivalent to the extrinsic curvature 
formulation. To see this, note that by virtue of (|2.56() . Eqs. (|2.54() and H2.55|l become 



A'3 = M'^ + — ^ 
2N 



h{V - W) 



and 



,27V 

With the identifications 



[L(y-w)y 



— —u 

2N 



(2.57) 



(2.58) 



(2.59) 



Eqs. H2.57|l and H2.58|l are identical to Eqs. (|2.24|) and (|2.25|) of the conformal thin 
sandwich formalism. The Lagrangian picture agrees completely with the Hamilton- 
ian picture. A third reason for (|2.56|) is given next. 
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2.3.3. Stationary spacetimes have A^^j. = 

Consider a stationary solution of Einstein's equations with timelike Killing vector I. 
Given a spacelike hypersurface S, there is a preferred gauge so that the time- vector 
of an evolution coincides with Z, namely N = —n • /, = ± /, where n is the unit 
normal to S, and _L is the projection operator into E. With this choice of lapse and 
shift, gij and K'^^ will be time-independent. Using dtgij = in Eq. H2.2|l and taking 
the tracefree part yields 



a weighted transverse traceless decomposition with = 0. Thus, with the ap- 
propriate weight factor a = 2A^, the extrinsic curvature has no transverse traceless 
piece for any spacelike slice in any spacetime with timelike Killing vector (A similar 
argument is applicable in the ergosphere of a Kerr black hole; however, one must 
be more careful with the choice of E relative to 

This is an important result. One generally identifies the transverse traceless piece 
of the extrinsic curvature with radiative degrees of freedom. Stationary spacetimes 
do not radiate, and therefore should indeed vanish. In contrast, a transverse- 
traceless decomposition of A^^ without the weight-factor will in general lead to a 
nonzero transverse traceless piece so that such a decomposition is incompatible with 
the identification of Al^j, with "gravitational radiation." 

3. Binary black hole initial data 

Up to about five years ago, many assumptions were necessary to simplify the initial 
value problem sufficiently to make it tractable for the computational methods of 
that time. The major assumptions were (see |12| for a review) 

(1) Maximal slicing, K = 0. 

(2) Conformal flatness gij = fij, where fij represents the Euclidean metric. 

(3) Use of the analytical Bowen-York |S] extrinsic curvature to solve the momentum 
constraint. 

Under these assumptions, only a single quasi-linear (flat-space) Laplace-equation 
must to be solved numerically for the conformal factor. This was done, e.g., with 
inversion symmetry boundary conditions |11I14| or with the puncture method 0. 
Since then, the conformal method (as presented in Sec. was completed, and 
several numerical codes were developed that are capable of solving the coupled 
constraint equations |31I23I25I38I46I47| . Three different numerical techniques have 
been used to discretize the elliptic problems, finite differences, spectral methods 
and finite elements. 

Finite differences 131146147) are familiar to almost all researchers, and are fairly 
straightforward to implement. However, the presence of different length scales 



(2.60) 



Ma,b ^d^R, 



(3.1) 
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where Ma,b represents the masses of the two black holes (labeled A and B), d their 
separation and R the distance to the outer boundary of the computational grid, 
restricts finite difference codes on uniform grids to a very coarse resolutions with 
errors of 10~^ to 10"'^. Imposing boundary conditions on spherical inner boundaries 
is also difficult in finite difference codes. Adaptive mesh-refinement may drastically 
improve the ability of finite-difference codes to handle the different length scales in 
Eq. 1)3.1(1 . In Refs. |17I41IH| the Hamiltonian constraint alone is solved, and work is 
in progress to extend adaptive mesh refinement to the coupled initial value equa- 
tions 42 . 

Spectral methods |23l38j (see |1()I7) for general introductions) approximate the 
solution with a truncated expansion in some basis functions, typically Chebyshev- 
polynomials or spherical harmonics. The solutions to the constraint equations are 
smooth, so that the accuracy improves exponential with the number of basis func- 
tions and much higher accuracies are achieved (~ 10~^ for binary black hole ini- 
tial data in |SH|)- Spectral methods are also more efficient, therefore permitting 
much larger parameter studies, allow more easily to impose boundary conditions on 
spherical boundaries (for expansions in spherical harmonics), and, in combination 
with domain-decomposition techniques, easily resolve the different length scales in 
Eq. (|3.1I) . The most accurate (for binary black holes) and versatile code seems to 
be the one developed by the author |2HI, which has been used to solve essentially 
all versions of the initial value problem |36l35ll5l?7| . 

Finite elements [25I4| cover the computational domain with very many small 
computational cells (typically tedrahedra), and expand the solution to low polyno- 
mial order (often just linear) in each cell. The method presented in Refs. [2514] is 
capable of solving the coupled constraint equations, and work toward using it for 
physically meaningful initial data is under way. 

With these new codes and with the final conformal method it is possible to 
move beyond the assumptions mentioned above. In particular, three separate issues 
have been pursued: Exploration of conformally non-flat three-geometries, which is 
motivated by the facts that the Kerr spacetime does not admit conformally flat 
slices |33l21l2y] . and that a binary compact object is not conformally flat at sec- 
ond post-Newtonian order 02] • Replacement of the Bowen- York extrinsic curvature, 
which does not exactly reproduce a stationary spinning or boosted black hole, and 
which may be responsible for unexpected behavior of sequences of circular orbits 
for spinning equal-mass black holes '39' or irrotational black holes in the test-mass 
limit 35 . And finally, physically motivated boundary conditions at inner excision 
surfaces surrounding the singularities of the black holes. 

The goal of these new investigations is the construction of astrophysically real- 
istic binary black hole data, i.e. configurations as they occur during the inspiral of 
two black holes in nature. Ideally, of course, such an initial data set should contain 
the outgoing gravitational wave signal of the preceding inspiral. However, at the 
current stage of sophistication, this is not taken into account, and the intermediate 
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goal is to construct initial data with as little "spurious," "unphysical" gravitational 
energy content as possible. The quotes used around the terms "spurious" and "un- 
physical" indicate part of the challenge: It is not even clear what these terms mean 
precisely on a single spacelike surface. 

Some of the research discussed in this section predates the conformal thin sand- 
wich formalism, or the final extrinsic curvature decomposition as discussed in Sec. [3 
These earlier papers typically use special cases of these more general frameworks, 
and we will discuss them from that perspective. 

3.1. Numerical solution of the extrinsic curvature decomposition 

As mentioned above, the Kerr spacetime does not admit conformally flat slices, so 
that any initial data built on conformal flatness will not be able to reproduce Kerr 
exactly. One approach to address this issue — proposed by Matzner et al. — 
superposes exact quantities for single (spinning or boosted) black holes to define 
the free data for the extrinsic curvature decomposition. Because the method is built 
on analytical single black hole solutions, it trivially works for single black holes. 
Specifically, Matzner et al. considered the Kerr-Schild form of the Kerr spacetime, 

^'^g^. = ^M- + 2i?y., (3.2) 

where rj^i, represents the Minkowski metric, if is a scalar function on spacetime, 
which decays as 1/r at large radn, and is null with respect to both the full metric 
and the Minkowski metric (the concrete expressions can be found, e.g., in |S21)- The 
Kerr-Schild form is preserved under a Lorentz-transformation, and the metric for 
a boosted black hole can be obtained by a suitable Lorentz-transformation on l^. 
In order to construct initial data for a spacetime containing two black holes at 
coordinate locations ca,b with masses Ma.b^ velocities va.b and spins Mao, a and 
MbSb, one constructs the Kerr-Schild form Eq. H3.2|l for each of the two black holes 
separately. The free data for the extrinsic curvature decomposition is then taken as 
the superposition 

ffy = % + 2H^lflf + 2i/^/f /f , (3.3) 
K^Ka + Kb, (3.4) 

M'' = + - l^^k (Ka + Kb)'^ ~g'^\ (3.5) 

Equation (|3.3|l is a natural generalization to two black holes of Eq. %\.2\ . Equa- 
tion l|3.5|) is somewhat complicated because M*-' must be tracefree with respect to 
the superposed metric gij, cf. Eq. (|2.51|l . In Ref. PJj the eUiptic equations of the 
extrinsic curvature decomposition were solved for the case d = 1 (i.e. one of the old 
variants of the extrinsic curvature decomposition). 

In the limit of large separation of the black holes, one obtains two widely sepa- 
rated Kerr-Schild metrics, each of which, by construction will represent a (possibly) 
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boosted and spinning black hole exactly, so that the superposed Kerr-Schild ap- 
proach certainly is advantageous for widely separated binaries. However, the most 
interesting initial data is for a binary with separation close to the innermost stable 
circular orbit, with black hole separations of only a few Schwarzschild radii. In that 
case, close to the horizon of, say, hole A, the correction 2H^lfl^ due to hole B 
will not be small; for example, at the intersection of the line connecting the cen- 
ters ca and cb with the horizon of hole A, one finds for separation lOM and for 
non-spinning black holes. 

Given the nonlincaritics in Einsteins equations, it is not clear how the final initial 
data set will be influenced by this large perturbation. Since it is expected that only 
a few percent of the energy will be emitted in gravitational waves during the inspiral 
and merger of a binary black hole, it is necessary to control the energy-content of the 
initial data set to at least one percent of the energy (and preferably much better). In 
view of the ratio (|3.6(l . it seems unlikely that this is the case. Indeed, the author [21] 
has examined the proposal of superposed Kerr-Schild quantities extensively with a 
spectral code for the construction of two black holes at rest. Using superposed 
Kerr-Schild quantities within the conformal thin sandwich formalism and within 
the extrinsic curvature formalism resulted in variations of the ADM-energy by up 
to several per cent. It was also found that the resulting initial data sets depend 
sensitively on the choice of extrinsic curvature: Removal of the longitudinal part of 
M^^ by a procedure analogous to Eqs. (|2.51l) - (|2.54|) before solving the constraint 
equations changes the ADM-energy by several percent. 

We also note that superposed Kerr-Schild metrics singles out arbitrarily a spe- 
cific slicing: Any slicing of the Kerr spacetime will give rise to three-metric and 
extrinsic curvatures which can be superposed similarly to Eqs. (|3.3|I - H3.5|I . Further- 
more, the superposition depends on the spatial coordinate system used to represent 
the single black hole quantities. For non-spinning, unboosted black holes, the spatial 
part of the Kerr-Schild metric is spherically symmetric, and can be made confor- 
mally flat by a radial coordinate transformation After this spatial coordinate 
transformation, the equivalent of Eq. (|3.3(l is superposition of two conformally flat 
metrics, which implies that g.ij should be chosen to be conformally flat |35I15| . While 
superposing single black hole quantities is certainly an interesting route, more ex- 
haustive investigations into their properties are necessary. 

Conformal flatness is also questionable because a binary compact object is not 
conformally flat at second post-Newtonian order Tichy et al. 0^1 address this 
issue by using post-Newtonian results to set the free data for the extrinsic curvature 
decomposition: The post-Newtonian expansion of the spatial metric can be written 
as 123 

5™ = ^lN%+/*r' (3-7) 
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where i/ipN and /i^-^ are given as expansions in the velocity v/c. The spatial metric 
Eq. (|3.7(l and the corresponding extrinsic curvature are known for general motion of 
the masses, i.e. in particular for a binary in circular orbits. The extrinsic curvature 
can be similarly expressed as a series in w/c, the leading order term being the 
Bowen-York momentum. One now bases the conformal metric on Eq. 1)3.7(1 . and 
similarly for the extrinsic curvature. This approach is the only one so far, which, in 
principle, can account for genuine binary contributions to the free data. 

Tichy et al. j46j encountered two important issues during implementation of this 
approach. The first problem is that close to the point-masses, the post-Newtonian 
expansion breaks down, that is, different terms in the expansion grow with dif- 
ferent inverse powers of distance to the point mass. Terms of higher order in v/c 
diverge faster than the lower order terms, so that the choice which terms one retains 
(e.g., consistent in the order v/c) influences the obtained initial data sets strongly. 
This may reflect a fundamental problem of using post-Newtonian expansions: They 
are least accurate close to the point-masses. The second problem encountered by 
Ref. is related to the fact that tpp-^ has a singularity at the point-masses. In 
order to obtain a finite conformal metric, the authors decided to conformally scale 



where the leading order behavior of fl close to the singularities is identical to that 
of V'PN However, the concrete choice for f2 influences the resulting initial data sets 
in a significant way; for example, the simplest choice — ^p-^ leads to binary black 
hole initial data sets with increasing energy as the separation between the black 
holes is reduced. Tichy et al. do not employ the extrinsic curvature decomposition 
presented in Sec. 12.31 but rather the old version without the weight-function a (their 
work predates discovery of the weight-function in Ref. 001 )• This old decomposition 
is not invariant to conformal transformations of the free data, while the new one 
is [cf. Eq. H2.50|l ]. so that ambiguities related to the choice of ft arise only in the 
old decomposition. It would be interesting to repeat Tichy's work with the new 
decomposition. 

3.2. Numerical solutions of the conformal thin sandwich equations 

The conformal thin sandwich formalism has two advantages over the extrinsic cur- 
vature decomposition. First, it replaces the free data corresponding to the trace- 
free extrinsic curvature M*-' and the weight function a by freely specifiable time- 
derivatives, dtQij = Uij and dtK. Time-derivatives allow for easier physical inter- 
pretation of the initial data set under construction, and as we will see below, in 
the most interesting case, there is a natural choice for these time-derivatives, thus 
avoiding ambiguities related to the choice of, e.g.. A/'-' in Eq. H3.5|l . As second 
advantage, solution of the conformal thin sandwich equations results directly in a 
preferred gauge choice (A^, /?*), which can be used in evolutions of the initial data 




n-4 




(3.8) 
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set (at least during the initial stages of the evolution). 

The orbits of binary compact objects of similar masses are expected to circularize 
before late inspiral, so that shortly before merger, the compact objects will move 
in circular orbits about each other. Therefore, initial data with two black holes in 
such a quasi-circular orbit is of particular interest. In the co- rotating coordinate 
system such a configuration will appear time-independent (up to corrections due to 
radiation reaction which are neglected so far), and so the free data corresponding 
to time-derivatives is simply set to zero. This motivates use of the conformal thin 
sandwich formalism with the free data Eq. H2.3()|l . fully half of which consists of 
time-derivatives: 

u„ = dtK = 0, (3.9) 

The notion of time-independence in the corotating frame is basically equivalent 
to the assumption of an approximate helical Killing vector |22], which close to 
infinity takes the form ^hkv — d/dto + flod/d(j>o, where d/dto and d/d(f>Q are, 
respectively, the asymptotic timelike and rotational Killing vectors. This follows 
simply from the fact that, in the corotating frame, the time- vector of the evolution 
coincides with the approximate helical Killing vector. Assuming asymptotic flatness, 
and requiring that the time-vector of a subsequent evolution coincides with ^hkv 
close to infinity leads in boundary conditions at the outer boundary. 



V'^1, r~~>oo, (3.10) 

P^^^o{^)\ .-oo, (3.11) 

iV^l, r->cx). (3.12) 



It now remains to choose gij and K as well as inner boundary conditions at 
excision spheres around the singularities of the black holes^. 

The first proposal for such boundary conditions was made by Gourgoulhon et 
al. |22|^. and numerical solutions were obtained with a spectral code in a companion 
paper Refs. |22I23| assumed conformal flatness and maximal slicing to simplify 
the problem. Part of the necessary boundary conditions at inner boundaries were 
derived from the demand that the obtained initial data be inversion symmetric 
across the throats of the horizons, resulting in a Robin boundary condition on the 
conformal factor, 

s'V.V + — = 0, on 5 (3.13) 
2r 

''The puncture method jSJ avoids inner boundaries in the computational domain, but it cannot be 
generaUzed to the construction of quasi-equiUbrium initial data within the conformal thin sandwich 
formalism 24 . 

''Gourgoulhon et al. scale the tracefree extrinsic curvature differently from Eq. 12.241 . so that the 
elliptic equations differ in lower order terms and some of the statements made in Sec. l2.2l do not 
apply. 
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where r denotes the coordinate radius of each excised sphere, and represents the 
outward pointing (i.e. toward infinity) normal to S, tangent to the hypersurface and 
normahzed such that gijS^s^ = 1. Inversion symmetry also implies that the lapse 
function must vanish on the throat, 

N = Q, on 5, (3.14) 

which constitutes a boundary condition on the fifth initial value equation, Eq. (|2.31|l . 
Boundary conditions on the shift, finally, were obtained by requiring that the time- 
vector i'* = Nn^^ + be tangent to the generators of the horizon. In particular, t^^ 
must be nuU, tH"'^'^'' g^^ = 0, so that from Eq. Ipnijl it follows 

(3' = 0, on S. (3.15) 

Since the lapse-function vanishes on S, inspection of Eqs. H2.23|l and 1)3. 9|l reveals 
that the extrinsic curvature will be finite on S only if (L/3)'^ vanishes. However, 
this is not the case, so that Refs. j22|23j must resort to a regularization procedure 
which causes the resulting initial data to violate the constraints on some small 
level. The authors continue and construct a sequence of quasi-circular orbits for 
corotating black holes, and find that the location of the innermost stable circular 
orbit (ISCO) is close to post-Newtonian predictions. 

A more general and more sophisticated approach to inner boundary conditions 
is due to Cook & Pfeiffer [13135115] . The idea is to use the physical concept of 
black holes in quasi-equilibrium to derive as many boundary conditions as possible. 
As we will see, this approach can be used in combination with any choice for gij 
and K. Two topological spherical regions with boundary S are excised from the 
computational domain. Denote the outward-pointing nuU-geodesics tangent to S 
by fc^, and their expansion by 6. Then the following demands are made: 

(1) The surfaces S must be apparent horizons'^, i.e. 9\g = 0. 

(2) The shear of A;'' must vanish on S. 

(3) When the initial data are evolved with the lapse and shift obtained during the 
solution of the conformal thin sandwich equations, then the coordinate surfaces 
of the apparent horizons are to remain stationary initially. 

The first of these demands simply localizes the apparent horizon in coordinate space, 
so that it is known where to apply the other conditions. The second requirement 
is based on quasi-equilibrium; it implies, using Raychaudhuri's equation for null- 
congruences, that 

CkOls^O. (3.16) 

That is, initially, the apparent horizon will evolve along and because the ex- 
pansion of vanishes, the apparent horizon area will remain constant initially. 
The third demand is a gauge choice ensuring that the coordinates are adapted to 



'^More precisely, marginally outer-trapped surfaces. 
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the physical situation. Reexpressing demands (1) and (3) in the variables of the 
conformal thin sandwich equations yields the following conditions: 

s^Vfcln^ = -2h''V^S, + ^^^K - ^g,5J(L/3)^^ on S, (3.17) 

= '4j^NS' + on 5, (3.18) 

where s*^ is defined below Eq. 1)3.13(1 . hij — cjij —SiSj denotes the induced conformal 
metric on S and is the projection of into S (i.e. = 0). Equation 1(3.17(1 
represents a nonlinear Robin-like boundary condition on -0, whereas Eq. ((3.18(1 is 
essentially a Dirichlet condition on the shift. The remaining demand of vanishing 
shear of fc^, finally, implies ^] that the tangential component of the shift /3*| must 
be a conformal Killing vector of the two-dimensional surface S, 

D^'pf - fh'WkPl =0 on 5, (3.19) 

where Di denotes the covariant derivative within S, which is induced by the con- 
formal metric. Cook & Pfeiffer point out that solutions to Eq. ((3.19(1 can be found 
before the conformal thin sandwich equations are solved: Each connected compo- 
nent of {S,hij) is topologically and therefore conformally equivalent to a unit 
2-sphere embedded in three-dimensional Euclidean space. For such a Euclidean 2- 
sphere, Euclidean rotations around the center represent Killing vectors, which are 
conformal Killing vectors of any conformally related manifold, including (iS, /ly). 
Therefore, these rotations will satisfy Eq. 1(3.19(1 . The freedom to specify arbitrary 
rotations can be used to construct black holes with arbitrary rotational state. In- 
terestingly, the lapse boundary condition is not determined by quasi-equilibrium 
considerations, but is rather part of the temporal gauge choice. 

The quasi-equilibrium boundary conditions at the inner excision regions allow for 
arbitrary specification of conformal metric, mean curvature, and shape of the exci- 
sion regions. Spectral numerical solutions of the conformal thin sandwich equations 
(with lapse equation) using the quasi-equilibrium boundary conditions Eqs. ((3.17(1 - 
((3.19(1 were obtained in for the special case of a conformally flat metric and 
excision of exact coordinate spheres. Three lapse boundary conditions at the ex- 
cised regions (Dirichlet, von Neumann and Robin), combined with irrotational and 
corotating black holes, and combined with two choices for the mean curvature K 
were explored, for a total of 12 sequences of quasi-circular orbits. Furthermore, spin- 
ning and boosted single black holes were constructed under the approximation of 
conformal flatness, and it was shown that the errors introduced by the use of a flat 
conformal metric are fairly small. Nonetheless, the freedom in the choice of gij, K, S 
and the lapse boundary condition can be used to further fine tune the method. 

There is a close connection between the quasi-equilibrium inner boundary con- 
ditions and the isolated horizon framework developed by Ashtekar and coworkers 
(see, e.g. (2I18I3( '). Jaramillo et al. "25" work out this relation very clearly. Not 
surprisingly their results are closely related to findings in Ref. (15( , including an in- 
dependent argument why the lapse boundary condition is not determined by quasi- 
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equilibrium/isolated horizon considerations. Ref. |26| does not present numerical 
results, and it is not immediately clear how to translate the conditions in this paper 
into a form usable in numerical simulations. Yo el. al. j47| recently solved the five 
coupled conformal thin sandwich equations on a Kerr-Schild background, however, 
with much simpler inner boundary conditions, and with much less accuracy, owing 
to the employed finite-difference code with uniform grid spacing. 

As a final, somewhat unrelated application of the conformal thin sandwich 
equations, we mention construction of initial data with superposed gravitational 
waves [HZj (see also j44l5) V The idea is very simple: Given a slice through a station- 
ary background spacetime with induced metric 5? and mean curvature , as well 
as a linearized gravitational wave hij (either on the background, or on Minkowski 
space), set the free data for the conformal thin sandwich equations as 

~g,,^gl+Ah,,, (3.20) 

u,^j ^ Adth^j, (3.21) 

K = (3.22) 

dtK = 0, (3.23) 

where A is the adjustable amplitude of the perturbation. In Rcf. ^37" the conformal 
thin sandwich equations were solved with these free data (and Dirichlct boundary 
conditions). Solutions were obtained for perturbations of Minkowski space and of 
a Schwarzschild black hole with very large amplitudes. These perturbed black hole 
initial data are used to test constraint preserving boundary conditions for evolutions 
of Einstein's equations |28| . 

4. Conclusion 

We summarized the conformal method to solve the constraints of general relativity, 
both in its Lagrangian and Hamiltonian viewpoints. By virtue of the recently dis- 
covered weighted transverse-traceless decomposition in the Hamiltonian viewpoint 
(the extrinsic curvature formulation), it is shown to be completely equivalent to the 
Lagrangian viewpoint (the conformal thin sandwich formulation). Both pictures are 
invariant to conformal transformations of the free data, and in both pictures, the 
decomposition of the tracefree extrinsic curvature commutes with the conformal 
scalings. 

Subsequently, we summarized recent numerical work on solving the coupled con- 
straint equations, in either viewpoint, concentrating on quasi-equilibrium solutions 
using the conformal thin sandwich equations, as well as superposed Kerr-Schild 
data and initial data based on post-Newtonian results using the extrinsic curvature 
formulation. 
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